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A continuous-time quantum walk is investigated on complex networks with the characteristic property of 
community structure, which is shared by most real-world networks. Motivated by the prospect of viable quantum 
networks, I focus on the effects of network instabilities in the form of broken links, and examine the response 
of the quantum walk to such failures. It is shown that the reconfiguration of the quantum walk is determined by 
the community structure of the network. In this context, quantum walks based on the adjacency and Laplacian 
matrices of the network are compared, and their responses to link failures is analyzed. 
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I. INTRODUCTION 

Networks are ubiquitous in both nature and society. They 
are routinely used to simulate a wealth of phenomena in the 
physical and biological sciences, as well as in sociology, fi- 
nance, information and communication technologies Jl], 0] ■ 
In the vast majority of such applications the employed net- 
works are inherently complex, by which we mean that there 
are strong fluctuations in their structural characteristics. This 
structural disorder is, in fact, a new type of disorder that can 
lead to cooperative behavior which goes beyond the one en- 
countered in traditional condensed matter physics |3]. 

Quantum networks have become a viable prospect in the 
area of quantum information processing, with potential ap- 
plications ranging from teleportation to cryptography 101 . In 
view of their potential use in the foreseeable future, it is 
clearly beneficial to determine the role of structural complex- 
ity in the dynamics of quantum networks. 

A small step in this direction is taken in the present work 
by focusing on a characteristic property of complex networks, 
which is typically referred to as community structure isj]. Intu- 
itively, a community is a cluster of nodes (vertices) in a com- 
plex network (graph), which is connected more densely on the 
inside than it is connected with the outside. In other words, 
there are more intracommunity links (edges) within the com- 
munity than there are intercommunity links between that par- 
ticular community and other communities in the network. 

As a straightforward illustration we shall examine a social 
network known as Zachary's karate club (KC) depicted 
in Fig. [T] The specific network has been studied extensively in 
the field of community detection i^. The main results are pre- 
sented here in relation to the KC network, but they are valid in 
general and apply equally well to other networks of increasing 
size and complexity. In particular, the results have been cor- 
roborated by calculations on the bottlenose dolphins network 
with iV = 62 nodes |@] and benchmark artificial networks of 
various sizes N G [40, 500] with heterogeneous community 
structure ifioll , studied for the altogether different purposes of 
community detection . 

In this setting, I examine a continuous-time quantum walk 
(CTQW) and its dynamical response to structural instabilities 
of complex networks. CTQWs have been studied well in dif- 
ferent contexts llTlllT2ll , including quantum search algorithms 
and quantum communication with spin system dynamics ifisll . 



CTQWs on statistical models, such as small-world networks, 
have also been studied lfT2ll . In the context of spin lattice 
dynamics ifTill and modified quantum walks, it was recently 
shown that quan tum walks can detect structural faults in reg- 
ular graphs Hslfiell. The aim of the present work is differ- 
ent, however, namely to examine the behavior of CTQWs on 
real-world networks and assess their behavior following a link 
failure (fault with the connections of the network). 




FIG. 1: (Color online) Community structure in the karate club (KC) 
network @] with A'^ = 34 nodes. The two main communities, cen- 
tered around nodes 1 and 34, are indicated by squares and circles, 
respectively. Colors correspond to the various possible communities 
in the network, including sub-communities. Reprinted from Ref. 0] 
((5)2004, lOP Publishing and SISSA). 

The rest of the paper is organized as follows. Sec. II in- 
troduces the model and defines the essential quantities to be 
used later on. Sec. Ill solves the model on the KC network 
and on larger systems, and establishes the main result, which 
can be quantified using the iwde affinity function, introduced 
here precisely for this purpose. Sec. IV extends the anal- 
ysis by comparing the behavior of the system with a differ- 
ent type of quantum walk, in order to assess the robustness 
in each case (in particular, CTQWs based on the adjacency 
and Laplacian matrices of a given network are compared and 
contrasted). Sec. V concludes with a summary of results, 
comments on experimental implementation, and potential ap- 
plications in quantum information science. 
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II. MODEL AND DEFINITIONS 

A complex network G{V, E), composed of = |V^| ver- 
tices and K ~ \E\ edges, can be described by an adjacency 
matrix A{Q), given by 



1 if (z,i)€i?(5), 



A,, 



Aij = otherwise 



for a network that is unweighted, undirected (Aij = Aji) and 
without loops (Ajj = 0). There are various other matrices 
that can be associated with a given network ITtIi . such as the 
Laplacian matrix 



L,. 



D, 



A,, 



where Dij = 6i,jdi is a diagonal matrix determined by the 
degree of each node j, that is, 

dj = ^ Aij . 

i 

With every node j = {1,2,..., N} of the network Q, we 
associate a basis state |j) in an iV-dimensional vector space. 
The basis states are orthonormal and a standard representation 
can be adopted, such as. 



|1> 



Any other state |?/') can then be written as a linear combina- 
tion, IV') = J2j CjIj), where Cj = 

At i = the initial state of the network is |*(0)). At later 
times, the evolution of the CTQW is given by 
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w 








[ij 



=exp(-iylt)|*(0)). 



(1) 



The evolution operator depends on the adjacency matrix A of 
the network. In the literature, by contrast, it is more common 
to use the Laplacian matrix instead lfTll[l2ll . Therefore, in the 
penultimate section of this paper, the two CTQWs are com- 
pared and contrasted, that is, in addition to the dynamics ob- 
tained from Eq. ([T]i we also examine the evolution according 
to the equation 



|vI'(t))L=exp(jLt)|vI'(0)). 



(2) 



Note that the subscript L on the state will be used to distin- 
guish it from the state of Eq. In this second case the evolu- 
tion operator can be decomposed into exp (iDt) exp {—iAt). 
Clearly, when the network is regular and each node has the 
same number of links to other nodes, i.e., dj = d for every j, 
the first term becomes a multiple of the unit matrix and there- 
fore the CTQWs given by Eqs. ([T]i and (|2]i are identical, up 
to an overall phase factor. However, in the case of complex 
(highly irregular) networks the two evolutions are different. 



The probability of finding the quantum walker on a node j 
at time t is 



P,(t)^|0-|vl,(t))|2 
and we have ^ Pjit) = 1- The time-averaged quantity 



Pjit)dt 



(3) 



(4) 



gives the mean probability of finding the walker on node j. In 
what follows, Pj plays a crucial role and it is referred to as the 
population of a node j. 

Finally, unless otherwise stated, we consider that there is an 
equal a priori probability to find the quantum walker on any 
node j at the start of evolution (at t — 0). This leads to an 
initial state of the form 



(5) 



In terms of a physical model, the adjacency matrix A is 
in fact the effective Hamiltonian of an XF-interaction spin 
model restricted to the single-excitation subspace. This model 
has been studied extensively in the field of quantum commu- 
nication with spin chains and lattices lUslQ. Therefore, in 
this context, the walker is a quantum excitation (i.e., a quasi- 
particle) diffusing in the network according to Eq. ^ and 
so, in the following, the term excitation is sometimes used to 
describe the quantum walker. 



III. MAIN RESULTS 

A. Populations Concentrate on Central Nodes 

One of the most basic structural characteristics of a com- 
plex network is the centrality of its nodes. A straightforward 
measure is the degree centrality. 



N - 1 



For Cj = the node is isolated; and for Cj = 1 the node is 
connected with every other node in the network. 

In order to probe the centrality of nodes using the formal- 
ism of CTQWs, we begin by assuming that there is equal a 
priori probability po = 1/iV of measuring the excitation on 
any given node. Therefore the initial state of the system is 
the state |^'(0)) of Eq. Q. Consequently, the evolved state 
l'I'(t)) of Eq. ([TJ is obtained numerically for a period of time 
t G [0, T], in time-steps St <C T, and finally the populations 
are calculated as prescribed by Eq. 

By comparing the population of each node with its degree 
centrality we see that populations generically tend to flow into 
highly-connected nodes. This is expected in the sense that, af- 
ter time t ^ N, the average probability of finding the walker 
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on a node j should increase in line with its centrality Cj: the 
more links that are incident on the node, the higher the proba- 
bility of the walker visiting that node. 

For the KC network in particular, the result is presented in 
Fig. 12] where both the degree centralities Cj and populations 
Pj are shown. It is clear that the population distribution is cor- 
related with the degree centrality of the nodes. The result has 
been corroborated by performing the same numerical calcula- 
tion on other real and artificial ifToll networks. The pop- 
ulation outcomes (average probabilities) remain stable with 
increasing integration times T, as long as T ~ cN, where c is 
a positive constant (typically c ~ 10). 




1 5 10 15 20 25 30 34 
Node 



FIG. 2: (Color online) Degree centrality Cj (line of squares) and 
population Pj (line of circles) for each node j — {1, 2, . . . , 34} 
in the KC network of Fig. [T] Parameters used for the numerical 
simulation: T = IOOtt (dimensionless units) and St = 10~"^r. 



B. Flow of Populations after a Link Failure 

We are now in a position to pose the central question of 
this work, namely. If a link fails, how do populations flow? 
In other words, the primary focus is on the reconfiguration of 
average probabilities of the CTQW following the failure (i.e., 
removal) of an edge from the network. 

We saw previously that populations tend to flow out of pe- 
ripheral nodes and into central nodes. So it is reasonable to 
expect that population hubs are formed around central nodes 
and that these hubs are sustained by intracommunity edges, 
connecting members of the community in which the hub be- 
longs. As a result, the removal of an edge strictly from the 
interior of a community should weaken the hub and therefore 
the population of the corresponding community as a whole 
should decrease (while of course the population of the rest of 
the network should increase). Another way to put the same 
idea is that the more edges there are inside a community, the 
more time spent by the quantum walker (e.g., excitation) in 
that community, and therefore the higher the probability of 
finding the walker in that region of the network. So if a link 



fails, this probability is reduced and it becomes more likely to 
find the quantum walker elsewhere (hence the flow of popula- 
tion out of the community). 
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FIG. 3: (Color online) Population flows after an edge is removed. 
Typical examples are shown: (a) edge (1, 3) belongs to the commu- 
nity centered around node 1 (squares in Fig. [T](; (b) edge (19, 33) 
belongs to the community centered around 34 (circles in Fig. [TJ; and 
(c) edge (3, 9) is an intercommunity edge. The (blue) line of circles 
corresponds to the populations of Fig. |2]for an ideal network, while 
the (green and red) diamonds correspond to the re-calculated (square 
and circle) populations after the link failure. 

This expectation has been tested and indeed verified numer- 
ically in the KC and other [3 networks by removing 
edges and recalculating populations after the removal opera- 
tions. Typically, if the edge belongs to a community, popula- 
tions flow out of that particular community and into neighbor- 



4 



ing ones. But if the edge removed connects two communities, 
the populations of both hubs inside these communities are in- 
creased, while the population of nodes close to the failed link 
(removed edge) are reduced. 

Three illustrative examples are presented in Fig. [3] In Fig. 
|3f a) mb)] the failed link belongs to the community indicated 
by squares (circles) in the KC network of Fig. [T] after its re- 
moval the populations of the squares (circles) are reduced, as 
seen by the green (red) diamonds. By contrast, in Fig. [Sjc) 
the failed link lies in-between the communities of squares and 
circles; and after its removal the population hubs in both com- 
munities are increased. 

Therefore, the answer to the main question is that the fail- 
ure of a Unk entails that populations flow out of the commu- 
nity in which it belongs and into neighboring ones. If the link 
is in-between communities, then the populations of both com- 
munity hubs are increased. 

This cooperative behavior, on the community level, has 
been corroborated with extensive numerical testing on other 
real ^ and artificial ifioll networks. The results are not pre- 
sented, as they do not add to the main argument, but they are 
broadly similar irrespectively of the total size of the network. 
However, in the case of overlapping communities in tailor- 
made networks lIToll the direction of population flow can be 
more ambiguous, in some special cases, such as those with 
very high inhomogeneity of the community distribution size. 
However, such cases are not particularly relevant for the po- 
tential applications of our proposal, as these are outlined in 
the Conclusions. 



C. Node Affinity Function 



for instance, it is seen that node 1 has high affinity with nodes 
2 to 8 and 11 to 14 (among others), while it has low affinity 
with nodes {9, 10, 15, 16, 19, 21} and 23 to 34. In some other 
cases, such as large networks with many inhomogeneous com- 
munities ifioll . the function does not reflect community struc- 
ture quite as clearly, but it still captures the community-based 
response to link failures, on average. However, the similarity 
of dynamical response of different nodes to structural instabil- 
ities, such as multiple link failures, is captured quite clearly. 




FIG. 4: (Color online) Node affinity mj of Eq. Q for the KC 
network of Fig. [T] Correspondence between shades (colors) and 
numerical values is shown in the vertical bar on the right. 



The previous results mean that the reconfiguration of pop- 
ulations following a link failure entails co-operative behavior 
on the community level. The populations inside a community 
flow in the same direction when an edge is removed. This 
motivates us to find a suitable measure of similarity between 
nodes, conveying the likelihood that two nodes will react in 
the same way in the event of a link failure. 

Therefore, we let the direction of flow on node j, after an 
edge fc G [1, K] has been removed, be given by 9j(k) = ±1, 
where +1 (—1) denotes population flowing into (out of) the 
node. We then define 

°^y=Z^ K ' 

k=l 

which takes values between —1 and 1. Whenever a^j > 
{aij < 0) the two nodes are (are not) likely to react in a 
similar way, while a^j = 1 {ajj = —1) implies that they 
are definitely (not) likely to react similarly. Clearly, is 
a function capturing the likelihood that two nodes belong to 
the same community and hence respond to link failures in the 
same way, so we call it node affinity. 

The node affinity function for the KC network is presented 
in Fig. |4] In this case, the function reflects the community 
structure very accurately indeed. Looking at the first column. 



IV. COMPARISON WITH LAPLACIAN-BASED CTQW 

We now turn our attention to the CTQW determined by the 
Laplacian matrix of the network, as prescribed by Eq. (O. The 
Laplacian matrix is also known as the connectivity matrix by 
some authors lfl2ll and the corresponding CTQW appears in 
various areas of physics, chemistry and biology. 

The spectrum of the Laplacian matrix can be derived from 
that of the adjacency matrix only for regular graphs. In the 
case of complex networks, such as those considered here, 
the spectra of the two matrices are different. In particular, 
the eigenvalue spectrum of the Laplacian matrix for a single- 
component network (i.e., one without isolated regions) is of 
the form 

ll<l2<--- <In, 

where the smallest eigenvalue is /i =0 since the matrix is 
positive semi-definite UtIi . Given that LI = 0, the eigen- 
vector corresponding to is the (normalized) vector 1 with 
all entries being equal to 1. This vector is, in fact, the state 
|5'(0)) given by Eq. (|5]i; in CTQW notation it is given by 
1 = J2j Therefore, the initial state |^(0)) will not evolve 
in time; we have |*(t))L = l*(0)) for all t. 
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Note that, conversely, if we consider regular graphs and we 
initiate the CTQW of Eq. ([T), which is governed by the adja- 
cency matrix, in the equiprobable state |^'(0)) of Eq. ([S]), then 
the system will not evolve. In this case also, the initial state 
is an eigenstate of the evolution operator (because the eigen- 
spectrum of the evolution operator based on A is the same as 
the one based on L for regular graphs). 

Consequently, in order to compare the two types of CTQW 
we need to start from a localized state, |^'(0)) = \j). Starting 
from this state, we calculate the adjacency-type CTQW of Eq. 
([T]i and the Laplacian-type CTQW of Eq. and then obtain 
the long-time average populations from and \'i>{t))-^, 

respectively, via Eq. Q. 

First, we see that in both cases the populations depend 
strongly on the initial state The diffusion process is quite 
different with this initial state, in the sense that much of the 
initial excitation remains localized in its starting point jlZ]. 

Second, in the case of the Laplacian-type CTQW, the pop- 
ulations derived from the long-time average of \'ii{t))-^ do not 
reflect the (degree) centrality of the nodes in the network. This 
is due to the strong dependence of the populations on the ini- 
tial state. To illustrate the point let us assume that we initi- 
ate the walk on a weakly-connected node with low centrality; 
then the long-time average probability of finding the walker 
on that node will still be high due to the initial conditions, 
even though the node in question has actually low centrality. 

Third, we find that the Laplacian-type CTQW does not re- 
spond to link failures (edge removal operations) on the level 
of community structure. In other words, the deletion of an 
edge from the network does not cause the populations inside 
the corresponding community to flow out of that community 
and into the rest of the network. Instead, we find that popula- 
tions flow in more complicated ways that do not correlate in 
general with the community structure of the network. 

V. CONCLUSIONS 

We have examined the dynamics of CTQWs on various 
complex networks [gl [l^l ■ When the evolution is governed 



by the adjacency matrix and the walk starts from an equiprob- 
able (delocalized) state, the population of each node reflects 
its centrality. When a link in the network fails (i.e., if an edge 
is removed) the populations reconfigure in a way that depends 
on the community structure of the network. 

In particular, if the failed link belongs to a community A, 
the populations inside A decrease while the populations out- 
side it increase. In other words, populations flow out of the 
community in which the link failure has taken place. By 
means of the node affinity function we have quantified the 
similarity of nodes in their response to such link failures. 

By contrast, for CTQWs based on the Laplacian matrix 
of the network, the dynamics is trivial if we start from the 
equiprobable state. In this case we need to initialize the walk 
in a localized state and, as a result, the correlation between 
node centrality and population is lost. The response to a link 
failure is not determined by the community structure of the 
network and so the node affinity function cannot detect simi- 
larities of response to link instabilities, among the nodes. 

There are various proposals for the experimental implemen- 
tation of CTQWs on (regular) networks, mostly with quantum 
optical methods [1 1,1 . Recent proposals for the simulation of 
exotic lattice systems with superconducting qubits could also 
be employed, especially for complex network structures, as 
they allow for arbitrary connectivity between sites 1 1 8(1 . 

The results of this work may have applications in quantum 
network design. The main idea would be that, alongside any 
other dedicated network function, quantum networks could be 
monitored in a continuous way by means of the CTQW pre- 
scribed by Eq. ([T]i. As long as the network is designed to 
have straightforward community structure (regular networks 
included) then the failure of a link, or the failure of neighbor- 
ing links within a finite area, would be detected by the popu- 
lation flows of the CTQW. 
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